Fast Implementation of Low Rank Approximation of a Sylvester Matrix
نویسندگان
چکیده
In [16], authors described an algorithm based on Structured Total Least Norm (STLN) for constructing a Sylvester matrix of given lower rank and obtaining the nearest perturbed polynomials with exact GCD of given degree. For their algorithm, the overall computation time depends on solving a sequence least squares (LS) problems. In this paper, a fast implementation for solving these LS problems is proposed. The increased efficiency is obtained by exploiting the low displacement rank of the involved coefficient matrices.
منابع مشابه
Fast Low Rank Approximation of a Sylvester Matrix by Structured Total Least Norm
The problem of approximating the greatest common divisor(GCD) for polynomials with inexact coefficients can be formulated as a low rank approximation problem with a Sylvester matrix. In this paper, we present an algorithm based on fast Structured Total Least Norm(STLN) for constructing a Sylvester matrix of given lower rank and obtaining the nearest perturbed polynomials with exact GCD of given...
متن کاملA Structured Rank - Revealing Method for Sylvester Matrix 1 )
We propose a displacement structure based rank-revealing algorithm for Sylvester matrix, then apply it to compute approximate greatest common division of two univariate polynomials with floating-point coefficients. This structured rank-revealing method is based on a stabilized version of the generalized Schur algorithm [8], and is a fast rank-revealing method in the sense that, all computations...
متن کاملExact Solutions in Structured Low-Rank Approximation
Structured low-rank approximation is the problem of minimizing a weighted Frobenius distance to a given matrix among all matrices of fixed rank in a linear space of matrices. We study the critical points of this optimization problem using algebraic geometry. A particular focus lies on Hankel matrices, Sylvester matrices and generic linear spaces.
متن کاملA Quadratically Convergent Algorithm for Structured Low-Rank Approximation
Structured Low-Rank Approximation is a problem arising in a wide range of applications in Numerical Analysis and Engineering Sciences. Given an input matrix M , the goal is to compute a matrix M ′ of given rank r in a linear or affine subspace E of matrices (usually encoding a specific structure) such that the Frobenius distance ‖M −M ′‖ is small. We propose a Newton-like iteration for solving ...
متن کاملComputing Approximate GCD of Univariate Polynomials by Structure Total Least Norm
The problem of approximating the greatest common divisor(GCD) for polynomials with inexact coefficients can be formulated as a low rank approximation problem with Sylvester matrix. This paper presents a method based on Structured Total Least Norm(STLN) for constructing the nearest Sylvester matrix of given lower rank. We present algorithms for computing the nearest GCD and a certified 2-GCD for...
متن کامل